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We show that low-angle grain boundaries (GB) in high-temperature superconductors exhibit 
intermediate Abrikosov vortices with Josephson cores, whose length I along GB is smaller that 
the London penetration depth, but larger than the coherence length. We found an exact solution 
for a periodic vortex structure moving along GB in a magnetic field H and calculated the flux 
flow resistivity Rf(H), and the nonlinear voltage-current characteristics. The predicted Rf(H) 
dependence describes well our experimental data on 7° unirradiated and irradiated Y BazCu^Oj 
bicrystals, from which the core size l(T), and the intrinsic depairing density Jb{T) on nanoscales 
of few GB dislocations were measured for the first time. The observed temperature dependence 
of JbiT) = Ji,o(l — T/T c ) 2 indicates a significant order parameter suppression in current channels 
between GB dislocation cores. 

PACS numbers: PACS numbers: 74.20. De, 74. 20. Hi, 74.60.-w 



Mechanisms of current transport through grain bound- 
aries (GB) in high-temperature superconductors (HTS) 
have attracted much attention, because a GB is a con- 
venient tool to probe the pairing symmetry of HTS by 
varying the misorientation angle & between the neighbor- 
ing crystallites jl], |^j. As $ increases, the spacing between 
the GB dislocations decreases, becoming comparable to 
the coherence length £ (T) at the angle $ ~ 4 — 6° . The 
exponential decrease of the GB critical current density 
■h = Jo ex p(— "iV^o) 01) makes GBs one of the princi- 
pal factors limiting critical currents of HTS || . Atomic 
structure of GBs revealed by high-resolution electron mi- 
croscopy have been used to determine local underdoped 
states of GB, defect- induced suppression of supercon- 
ducting properties at the nanoscale and controlled in- 
crease of Jb by overdoping of GB|^|, |j. Much progress 
has been made in understanding the microscopic factors 
controlling J,($) at zero magnetic field, but the behavior 
of vortices on GBs is known to much lesser extent. 

The extreme sensitivity of Jb('&) to the misorientation 
angle makes GB a unique tool to trace the fundamen- 
tal transition between Abrikosov (A) and Josephson (J) 
vortices Q in a magnetic field H above the lower critical 
field H c i. For $ <C ■do, vortices on a GB are A vortices 
with normal cores pinned by GB dislocations Q. For 
■& > i?oi thfi maximum vortex current density circulat- 
ing across the GB is limited to its intrinsic Jb($), much 
smaller then the bulk depairing current density Jd- Be- 
cause vortex currents must cross the GB which can only 
sustain Jb <C Jd, the normal core of an A vortex turns 
into a J core, whose length I ~ £Jd/Jb along the GB is 
greater then £, but smaller then the London penetration 



depth A, if J b > J<i/k, where k = A/£ ~ 10 2 f|. As <d 
increases, the core length 1(d) ~ ^Jd/Jbi^) increases, so 
the GB vortices evolve from A vortices for <C i)o to 
mixed Abrikosov vortices with Josephson cores (AJ vor- 
tices) at Jd/ 'k < Jb(i?) < Jd- The AJ vortices turn into 
J vortices at higher angles, ■& > ■dj ~ #o ln(«Jo/ Jd), for 
which l($) exceeds A. For $ = 5°, k = 100, and Jo = Jd, 
the AJ vortices determine the in-field behavior of GBs in 
the crucial region < ■& < i?j ~ 23° of the exponential 
drop of Jh(i?) (in a film of thickness d«A, the A J region 
■& < dj ~ $ m (A 2 /<i£) broadens even further). 

The A J structures have two length scales: the core size 
I > £ and the intervortex spacing a = (<f>o/ 'B) 1 / 2 . The 
larger core of AJ vortices leads to their weaker pinning 
along a GB, which thus becomes a channel for motion 
of AJ vortices between pinned A vortices in the grains 
H 0] (Fig. 1). The percolative motion of A J vortices 
along GBs gives rise to a linear region in the V — I 
characteristic of HTS polycrystals|| |[ |j| 10' 0- How- 
ever no present experimental techniques can probe the 
cores of GB vortices, because the lack of the normal core 
makes A J vortices "invisible" under STM, while neither 
the Lorentz microscopy nor magneto-optics have suffi- 
cient spatial resolution to distinguish A and AJ vortices. 
In this Letter we report a combined theoretical and ex- 
perimental analysis which enabled us to prove the ex- 
istence of AJ vortices in 7° Y Ba-iCu^O'j bicrystals and 
extract the core size l(T), and the intrinsic depairing cur- 
rent density Jb at the GB from transport measurements. 
The value Jb of a GB turns out to be much higher than 
its global critical current density J q b , which is limited 
by self-field and pinning effects g [n| [jj> Thc 
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FIG. 1: Current streamlines around A J vortices on a GB 
(dashed line) and the bulk A vortices in the grains, calculated 
from Eq. (5) for I = 0.2a. 



field region in which only a single AJ vortex chain moves 
along GB while the bulk A vortices remain pinned, can 
be considerably expanded by irradiation, as shown below. 

For H > H cl , both I = £J d /J b , and a = {fa/B) 1 / 2 
are smaller than A, thus the AJ vortices on low-angle 
GBs are described by a nonlocal equation for the phase 
difference 9{x,t) in the overdamped limit || 



1'(u)du 



I = C(j)o/16TT 2 X 2 J b 



sin 61 + /3, (1) 
T = 0o/27rci?J 6 , (2) 



where the overdot and the prime denote differentiation 
with respect to the time t and the coordinate x along GB, 
R is the quasiparticle resistance of GB per unit area, <fio 
is the flux quantum, c is the speed of light, /3 = J/Jb, 
and J(x) is the current density through GB induced by 
A vortices. Here ft = [3o + 5f3(x) is a sum of the constant 
transport current (3q and an oscillating component 8(3(x) 
due to the discreteness of the A vortex lattice. The term 
6f3(x) gives rise to a critical current J gb of the GB due 
to pinning of AJ vortices by A vortices in the grains |t], 
Q. Eqs. (Q) and (||) are independent of the pairing 
symmetry (which only affects J b ) and are valid for both 
bulk samples and thin films in a perpendicular field. 

We consider a rapidly moving AJ structure in the flux 
flow state, [3 3> /3 G! f° r which the pinning term 8j3(x) <C 
1 can be neglected O]. In this case Eq. (|l|) has the 
following exact solution that describes a stable periodic 
vortex structure moving with a constant velocity v: 

6 = 7T + 7 + 2 tan -1 [M tan k(x - vt) /2] , (3) 

s 2 = [y/(l - (3 2 + h) 2 + 4f3 2 h -l-h + /3 2 }/2h (4) 



Here s 



v/vq, vq = l/r, tan 7 
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FIG. 2: The V — J curves calculated from Eq. (6) for different 
magnetic fields h = H/Hq: 0.01(1), 0.05(2), 0.1(3), 0.5(4), 
10(5). Inset shows the field dependence of Rf(H). 



M = [1 + l/h(l + s 2 )] 1 ' 2 + [h(l + s 2 )]- 1 ' 2 , k = 2n/a, and 
a is the period of the AJ structure. For a/l — > oo, Eq. (|^) 
describes a moving chain of single A J vortices §|. Gen- 
erally, a(H) is different from the period of the A lattice, 
but for H ^> H c i, the spacing a = (cfro/H) 1 / 2 is fixed by 
the flux quantization condition. Eq. (^|) corresponds to 
the following field distribution H(x, y) produced by AJ 
vortices in the region \y\ < A: 



H = 



2ir\ 2 



Relnsin[a; — vt + i(\y\ + yo)]; 



(5) 



where sinhfcyo = \/h(l + s 2 ). Eq. (g) satisfies the 
Maxwell equation V 2 H = with the boundary condi- 
tion H 1 = (47r/c)[J b sin6» - hv0'/2eR — J] on GB, where 
6{x — vt) is given by Eq. (||). Fig. 1 shows the current 
streamlines calculated from Eq. (||). For y > 0, these 
streamlines coincide with those of a chain of moving fic- 
titious A vortices displaced by y — —yo away from GB. 

The mean voltage V on a GB is determined by the 
Faraday law, V — 4>ov/ca, which yields 



V = V 



(l-/3o 2 



hf 



1 



1/2 



(6) 



where V = RJ b /V2. The V - J curve shown in Fig. 2 
is similar to that obtained by molecular dynamic simu- 
lations of incommensurate vortex channels p^ j. In our 
case the nonlinearity of V(J) is due to the A J core ex- 
pansion as J increases ||. For J « Jj. the V-J curve is 
linear, V = RfJ, where Rf = Rwh/Jl + h) is the flux 
flow resistivity due to the viscous motion of AJ vortices. 
If H > Hd, then h = (2nl/a) 2 = H/H , and 



Rf = 



rVh 

VH + H a 



Hn = 



(2tt0 2 



(7) 
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FIG. 3: Rf(H) data points extracted from the slopes of 
the V(J) curves at 77K and different H for unirradiated 7° 
bicrystal. The solid curves are described by Eq. (7) with 
R — 4.05mQ and Ho = 0.14T. Inset shows V — J curves for 
0.17kOe, 0.51kOe, 0.85kOe, 1.7kOe, and 3.4kOe (from bot- 
tom to top, respectively). 

At H < H , Eq. (0) describes R F (H) for A J vor- 
tices, whose cores do not overlap. In this case Rf(H) is 
reminiscent of the ID Bardecn-Stcphcn formula, Rbs — 
R^/H/H C 2, except that in Eq. (Q) the core struc- 
ture is taken into account exactly. For H > Hq ~ 
(Jb/ Jd) 2 H c2 <C H c2 , the AJ cores overlap, and Eq. (Q) 
describes a crossover to a field-independent quasiparti- 
cle resistance of GB. This regime has no analogs for A 
vortices, whose normal cores overlap only at H c2 . The 
simplicity of Eq. (Q) enabled us to extract intrinsic GB 
properties from the measurements of Rp{H). 

We observed the AJ vortex behavior on 7° 
y ' Ba 2 CuzO-i , bicrystals with a sharp resistive transition 
AT < OAK at T c = 91K. Thin films of thickness 250 
nm were grown on [001]-oriented SrTiO^ bicrystals by 
pulsed laser deposition at 210 mTorr oxygen pressure 
and 810°C, and then annealed in oxygen at 830 Torr and 
520°C for 30 min. One 7° bicrystal was irradiated with 
lGeV Pb ions at a fluence corresponding to IT. Bridges 
25 [im, wide were patterned by Ar ion beam etching on 
a cooled sample mount to produce a four-point measure- 
ment geometry, as described in Ref. 0. The voltage 
probes were lOO^m apart, on cither sides of the GB. 
V — I curves were measured in a gas-flow cryostat in 
fields < H < 10 T. Both samples are in the transition 
region of $ between strongly coupled bicrystals that be- 
have in a manner similar to the single-crystal grains, and 
weakly coupled bicrystals that demonstrate the Joseph- 
son effect J13]] . The intragrain J c values at 77 K were 
0.1M A/ crnr and 0.27A/ A/cm 2 for the unirradiated and 
irradiated samples, respectively. 
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FIG. 4: Rf{H) extracted from the slopes of V(J) for dif- 
ferent T and H for the irradiated 7° bicrystal. The V — J 
curves have the extended linear regions similar to those for 
the unirradiated sample shown in Fig. 3. The solid curves 
are described by Eq. (7) with two fit parameters 7? and Ho- 



Fig. 3 shows V-I characteristics for the unirradiated 
sample at 77K and different magnetic fields. The V-I 
curves exhibit linear flux flow portions in a wide range 
of / above the depinning critical current I g b(H) which 
decreases with H p3| . In this work we focus on the 
flux flow region / ^> I g i,, where the dynamic resistance 
Rf(H) — dV/dl increases as \f~H at low H, but lev- 
els off at higher H. We found that Eq. (Q) describes 
the observed Rf(H) very well, thus the vortex cores 
on this GB overlap at H ~ Hq, well below H c2 . Be- 
cause the GB can sustain a finite supercurrent I g b even 
for H > Ho, the GB vortices lack normal cores. The 
good fit in Fig. 3 enabled us to extract R = 4.05mf2 
and H = 0.UT < H c2 by plotting 1/R 2 F (H) versus 
1/H, as shown in Fig. 2. Using Eq. we obtain 

that I = {fa/ Ho) 1 / 2 /2tt = 190A at 77K, thus vortices 
on this 7° GB are indeed AJ vortices with phase cores 
much greater than £(T) = ^/^/^-T/T c ~ 40l but 
smaller than A = - T/T c ~ 4000 A for £ - 15 A, 

A ~ 1500A, and T c = 91K. 

Rf(H) data for the irradiated bicrystal at different T 
are shown in Fig. 4. The good agreement between Eq. 
(0) and the observed Rf{H) enabled us to extract the 
temperature dependences of Ho and R shown in Fig. 5. 
While R(T) varies only weakly, the field Hq{T) exhibits 
a parabolic dependence, H (T) = H (0)(1 - T/T c f with 
H o (0) ~ 42T. As follows from Eq. (g), the fact that 
H (T) cx (T c - T) 2 implies l(T) oc (T c - T) -1 , which, in 
turn, indicates the SNS behavior of J b (T) cx (T c — T) 2 , 
if A(T) cx (T c - T)- 1 ' 2 in Eq. (§). From the data in 
Fig. 5, we can also obtain the absolute value of the 
intrinsic depairing current density J& averaged over the 
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FIG. 5: R(T) and Ho(T) extracted from the fit procedure 
represented in Fig. 4. The dashed line shows Ho = 42(f — 
T/T c ) 2 [Tesla]. 



Josephson core length I. To do so, we write Eq. (|7|) in 
the form J b = (3 J<j/4) 6nHo /H c2 , which express J b in 
terms of the measured parameters -Ho, H c2 = </>o/27r£ 2 , 
and J d = c0 o /12V37r 2 A 2 C. For H c2 (T) = H' c2 {T c - T), 
this yields J b ~ (3J d /4)[6 7 r J ffo(0)/T c ^ 2 ] 1 /2(i_ T / Tc )i/2 j 

whence J b (85K) ~ 0.3J d (85i4T) for F o (0) = 42T, and 
= 2T/i£'. Likewise we get J b (77K) ~ 0.23J d (77JC) 
for the unirradiated sample, (H (77K) = 0.14T). There- 
fore, our data indicate a significant suppression of the 
order parameter, even on the low-angle 7° GB, in agree- 
ment with the model of Ref . jl6| . 

For the observed H (T), the A J core length l(T) = 
(4> /H (T)) 1 / 2 /2ir ~ ll jl-T/T c)- 1 ^], ex ceeds the bulk 
coherence length £(T) = £ /y/l-T/T c at T c - T < 
T c , but remains smaller than A(T), except very close 
to T c . For instance, we obtain Z(80ii") ~ 91 A, while 
£(80*0 ~ 43A, and X(80K) = 4300l. The length 1{T) 
also exceeds the GB dislocation spacing w 32A, thus the 
moving AJ cores probe GB properties averaged over few 
current channels between dislocations. The core length 
~ £,Jd/Jb($) increases as ■& increases, becoming 
larger than A, if $ > i?o In K, in which case A J vortices 
turn into J vortices. Since the ratio J b /Jd decreases as 
T increases, higher-angle GBs can exhibit AJ vortices at 
low T and J vortices at T m T c , while lower angle GBs 
have A vortices at low T and A J vortices at T sw T c . 

The depinning current density J gb seen in Fig. 3 is 
due to interactions of AJ vortices with microstructural 
inhomogeneities along GB and pinned A vortices in the 
grains Q. Another essential mechanism occurs if the pe- 
riods of A J and A vortices do not coincide, resulting in 
misfit dislocations in the A J chain fll5[. These vortex dis- 
locations can strongly limit J gb because their intrinsic 
pinning by the A vortex lattice is exponentially weakjl5[. 
For H ^> H&, when the A J and A periods are close, the 
pinning of a few misfit dislocations in the AJ chain may 
be thus dominated by macroscopic T c and J b (x) varia- 



tions along GB caused by facet structures, strain fields, 
local nonstoichiometry, etc. JlTj]. In this Letter we focus 
on the flux flow state J> J gb , where pinning weakly 
affects Ftp. In this case measurements of Rf(H) reveal 
the essential physics of GB vortices whose moving AJ 
cores probe intrinsic properties of GBs at the nanoscale, 
regardless of particular pinning mechanisms which are es- 
sential at J J gb . Because the observed V — J curves are 
nearly linear above J gb (77K, IT) ~ 10 4 - 10 5 A/ cm 2 < 
J b ~ (0.2 - 0.3) J d ~ 1 - 10MA/cm 2 , the pinning region 
J ~ J gb is much smaller than the scale of Fig. 2. Thus, 
the nonlinearitry of V(J) due to the A J core expansion 
does not affect Rp measured at J < 3J gb so Eq. (0) can 
be used to fit the data. A similar approach was used 
to measure the flux flow resistivity of pinned A vortices 
driven by strong current pulses well above J c |fj"8[ . 

The case when only a single AJ vortex row moves along 
the GB, while the intragrain A vortices remain pinned 
corresponds to low fields, H < Hi. For H > Hi, the 
AJ vortices drag neighboring A vortices in the flux flow 
channel along GB § . The field Hi can be estimated from 
the condition that the pinning force / = (f>oAH/2a(H) 
of AJ vortices due to their magnetic interaction with 
A vortices equals the bulk pinning force (f>Qj c /c, where 
AH = <j) e~ 27:u / a /tt\ 2 is the amplitude of the oscillating 
part of the local field H(x) m B + AH cos(2ttx / a) pro- 
duced by A vortices along GB, and u is the spacing of 
the first A vortex row from GB[Q. Therefore, 

Hi~^J 2 c {Hi)/^AH 2 (8) 



Transitions from a single to a multiple row vortex mo- 
tion H result in a sharp upturn of the Rf (B, 77 K) 
curve at Hi ~ 4T in Fig. 4. Here Hi for our unirra- 
diated sample (J c = 0.1MA/cm 2 ) is 7.3 times smaller 
than Hi for the irradiated one (J c = 0.27 A/ A/cm 2 ). 
For X(77K) ~ 4000A, and Hi = AT, Eq. (|) yields 
u = ln(c 2 J ffi?!)o/47r 2 A 4 J 2 )/47r = 0.92a. 

In conclusion, vortices on low-angle grain boundaries 
in HTS are mixed Abrikosov- Josephson vortices. Exact 
solutions for a moving AJ vortex structure, the nonlinear 
V — J characteristic and the flux flow resistivity Rp(B) 
were obtained. From the measurements of Rf(B) on a 7° 
YBCO bicrystal, we extracted the length of the AJ core 
and the intrinsic depairing current density J b on GB. The 
analysis proposed in this work can be used for systematic 
studies of the effect of overdoping Q on the local sup- 
pression of the order parameter and current transport 
through nanoscale channels between GB dislocations. 
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